The Brain Gene Registry: a data snapshot.

Monogenic disorders account for a large proportion of population-attributable risk for neurodevelopmental disabilities. However, the data necessary to infer a causal relationship between a given genetic variant and a particular neurodevelopmental disorder is often lacking. Recognizing this scientific roadblock, 13 Intellectual and Developmental Disabilities Research Centers (IDDRCs) formed a consortium to create the Brain Gene Registry (BGR), a repository pairing clinical genetic data with phenotypic data from participants with variants in putative brain genes. Phenotypic profiles are assembled from the electronic health record (EHR) and a battery of remotely administered standardized assessments collectively referred to as the Rapid Neurobehavioral Assessment Protocol (RNAP), which include cognitive, neurologic, and neuropsychiatric assessments, as well as assessments for attention deficit hyperactivity disorder (ADHD) and autism spectrum disorder (ASD). Co-enrollment of BGR participants in the Clinical Genome Resource's (ClinGen's) GenomeConnect enables display of variant information in ClinVar. The BGR currently contains data on 479 participants who are 55% male, 6% Asian, 6% Black or African American, 76% white, and 12% Hispanic/Latine. Over 200 genes are represented in the BGR, with 12 or more participants harboring variants in each of these genes: CACNA1A, DNMT3A, SLC6A1, SETD5, and MYT1L. More than 30% of variants are de novo and 43% are classified as variants of uncertain significance (VUSs). Mean standard scores on cognitive or developmental screens are below average for the BGR cohort. EHR data reveal developmental delay as the earliest and most common diagnosis in this sample, followed by speech and language disorders, ASD, and ADHD. BGR data has already been used to accelerate gene-disease validity curation of 36 genes evaluated by ClinGen's BGR Intellectual Disability (ID)-Autism (ASD) Gene Curation Expert Panel. In summary, the BGR is a resource for use by stakeholders interested in advancing translational research for brain genes and continues to recruit participants with clinically reported variants to establish a rich and well-characterized national resource to promote research on neurodevelopmental disorders.

Clinical variants paired with phenotype: A rich resource for brain gene curation.

PURPOSE: Clinically ascertained variants are under-utilized in neurodevelopmental disorder research. We established the Brain Gene Registry (BGR) to coregister clinically identified variants in putative brain genes with participant phenotypes. Here, we report 179 genetic variants in the first 179 BGR registrants and analyze the proportion that were novel to ClinVar at the time of entry and those that were absent in other disease databases. METHODS: From 10 academically affiliated institutions, 179 individuals with 179 variants were enrolled into the BGR. Variants were cross-referenced for previous presence in ClinVar and for presence in 6 other genetic databases. RESULTS: Of 179 variants in 76 genes, 76 (42.5%) were novel to ClinVar, and 62 (34.6%) were absent from all databases analyzed. Of the 103 variants present in ClinVar, 37 (35.9%) were uncertain (ClinVar aggregate classification of variant of uncertain significance or conflicting classifications). For 5 variants, the aggregate ClinVar classification was inconsistent with the interpretation from the BGR site-provided classification. CONCLUSION: A significant proportion of clinical variants that are novel or uncertain are not shared, limiting the evidence base for new gene-disease relationships. Registration of paired clinical genetic test results with phenotype has the potential to advance knowledge of the relationships between genes and neurodevelopmental disorders.

Dominant negative variants in KIF5B cause osteogenesis imperfecta via down regulation of mTOR signaling.

BACKGROUND: Kinesin motor proteins transport intracellular cargo, including mRNA, proteins, and organelles. Pathogenic variants in kinesin-related genes have been implicated in neurodevelopmental disorders and skeletal dysplasias. We identified de novo, heterozygous variants in KIF5B, encoding a kinesin-1 subunit, in four individuals with osteogenesis imperfecta. The variants cluster within the highly conserved kinesin motor domain and are predicted to interfere with nucleotide binding, although the mechanistic consequences on cell signaling and function are unknown. METHODS: To understand the in vivo genetic mechanism of KIF5B variants, we modeled the p.Thr87Ile variant that was found in two patients in the C. elegans ortholog, unc-116, at the corresponding position (Thr90Ile) by CRISPR/Cas9 editing and performed functional analysis. Next, we studied the cellular and molecular consequences of the recurrent p.Thr87Ile variant by microscopy, RNA and protein analysis in NIH3T3 cells, primary human fibroblasts and bone biopsy. RESULTS: C. elegans heterozygous for the unc-116 Thr90Ile variant displayed abnormal body length and motility phenotypes that were suppressed by additional copies of the wild type allele, consistent with a dominant negative mechanism. Time-lapse imaging of GFP-tagged mitochondria showed defective mitochondria transport in unc-116 Thr90Ile neurons providing strong evidence for disrupted kinesin motor function. Microscopy studies in human cells showed dilated endoplasmic reticulum, multiple intracellular vacuoles, and abnormal distribution of the Golgi complex, supporting an intracellular trafficking defect. RNA sequencing, proteomic analysis, and bone immunohistochemistry demonstrated down regulation of the mTOR signaling pathway that was partially rescued with leucine supplementation in patient cells. CONCLUSION: We report dominant negative variants in the KIF5B kinesin motor domain in individuals with osteogenesis imperfecta. This study expands the spectrum of kinesin-related disorders and identifies dysregulated signaling targets for KIF5B in skeletal development.

Macrocephaly and developmental delay caused by missense variants in RAB5C.

Rab GTPases are important regulators of intracellular vesicular trafficking. RAB5C is a member of the Rab GTPase family that plays an important role in the endocytic pathway, membrane protein recycling and signaling. Here we report on 12 individuals with nine different heterozygous de novo variants in RAB5C. All but one patient with missense variants (n = 9) exhibited macrocephaly, combined with mild-to-moderate developmental delay. Patients with loss of function variants (n = 2) had an apparently more severe clinical phenotype with refractory epilepsy and intellectual disability but a normal head circumference. Four missense variants were investigated experimentally. In vitro biochemical studies revealed that all four variants were damaging, resulting in increased nucleotide exchange rate, attenuated responsivity to guanine exchange factors and heterogeneous effects on interactions with effector proteins. Studies in C. elegans confirmed that all four variants were damaging in vivo and showed defects in endocytic pathway function. The variant heterozygotes displayed phenotypes that were not observed in null heterozygotes, with two shown to be through a dominant negative mechanism. Expression of the human RAB5C variants in zebrafish embryos resulted in defective development, further underscoring the damaging effects of the RAB5C variants. Our combined bioinformatic, in vitro and in vivo experimental studies and clinical data support the association of RAB5C missense variants with a neurodevelopmental disorder characterized by macrocephaly and mild-to-moderate developmental delay through disruption of the endocytic pathway.

PSMC3 proteasome subunit variants are associated with neurodevelopmental delay and type I interferon production.

A critical step in preserving protein homeostasis is the recognition, binding, unfolding, and translocation of protein substrates by six AAA-ATPase proteasome subunits (ATPase-associated with various cellular activities) termed PSMC1-6, which are required for degradation of proteins by 26S proteasomes. Here, we identified 15 de novo missense variants in the PSMC3 gene encoding the AAA-ATPase proteasome subunit PSMC3/Rpt5 in 23 unrelated heterozygous patients with an autosomal dominant form of neurodevelopmental delay and intellectual disability. Expression of PSMC3 variants in mouse neuronal cultures led to altered dendrite development, and deletion of the PSMC3 fly ortholog Rpt5 impaired reversal learning capabilities in fruit flies. Structural modeling as well as proteomic and transcriptomic analyses of T cells derived from patients with PSMC3 variants implicated the PSMC3 variants in proteasome dysfunction through disruption of substrate translocation, induction of proteotoxic stress, and alterations in proteins controlling developmental and innate immune programs. The proteostatic perturbations in T cells from patients with PSMC3 variants correlated with a dysregulation in type I interferon (IFN) signaling in these T cells, which could be blocked by inhibition of the intracellular stress sensor protein kinase R (PKR). These results suggest that proteotoxic stress activated PKR in patient-derived T cells, resulting in a type I IFN response. The potential relationship among proteosome dysfunction, type I IFN production, and neurodevelopment suggests new directions in our understanding of pathogenesis in some neurodevelopmental disorders.

Mechanism of KMT5B haploinsufficiency in neurodevelopment in humans and mice.

Pathogenic variants in KMT5B, a lysine methyltransferase, are associated with global developmental delay, macrocephaly, autism, and congenital anomalies (OMIM# 617788). Given the relatively recent discovery of this disorder, it has not been fully characterized. Deep phenotyping of the largest (n = 43) patient cohort to date identified that hypotonia and congenital heart defects are prominent features that were previously not associated with this syndrome. Both missense variants and putative loss-of-function variants resulted in slow growth in patient-derived cell lines. KMT5B homozygous knockout mice were smaller in size than their wild-type littermates but did not have significantly smaller brains, suggesting relative macrocephaly, also noted as a prominent clinical feature. RNA sequencing of patient lymphoblasts and Kmt5b haploinsufficient mouse brains identified differentially expressed pathways associated with nervous system development and function including axon guidance signaling. Overall, we identified additional pathogenic variants and clinical features in KMT5B-related neurodevelopmental disorder and provide insights into the molecular mechanisms of the disorder using multiple model systems.

GABBR1 monoallelic de novo variants linked to neurodevelopmental delay and epilepsy.

GABAB receptors are obligatory heterodimers responsible for prolonged neuronal inhibition in the central nervous system. The two receptor subunits are encoded by GABBR1 and GABBR2. Variants in GABBR2 have been associated with a Rett-like phenotype (MIM: 617903), epileptic encephalopathy (MIM: 617904), and milder forms of developmental delay with absence epilepsy. To date, however, no phenotypes associated with pathogenic variants of GABBR1 have been established. Through GeneMatcher, we have ascertained four individuals who each have a monoallelic GABBR1 de novo non-synonymous variant; these individuals exhibit motor and/or language delay, ranging from mild to severe, and in one case, epilepsy. Further phenotypic features include varying degrees of intellectual disability, learning difficulties, autism, ADHD, ODD, sleep disorders, and muscular hypotonia. We functionally characterized the four de novo GABBR1 variants, p.Glu368Asp, p.Ala397Val, p.Ala535Thr, and p.Gly673Asp, in transfected HEK293 cells. GABA fails to efficiently activate the variant receptors, most likely leading to an increase in the excitation/inhibition balance in the central nervous system. Variant p.Gly673Asp in transmembrane domain 3 (TMD3) renders the receptor completely inactive, consistent with failure of the receptor to reach the cell surface. p.Glu368Asp is located near the orthosteric binding site and reduces GABA potency and efficacy at the receptor. GABA exhibits normal potency but decreased efficacy at the p.Ala397Val and p.Ala535Thr variants. Functional characterization of GABBR1-related variants provides a rationale for understanding the severity of disease phenotypes and points to possible therapeutic strategies.

Functional analysis of a novel de novo variant in PPP5C associated with microcephaly, seizures, and developmental delay.

We describe a proband evaluated through the Undiagnosed Diseases Network (UDN) who presented with microcephaly, developmental delay, and refractory epilepsy with a de novo p.Ala47Thr missense variant in the protein phosphatase gene, PPP5C. This gene has not previously been associated with a Mendelian disease, and based on the population database, gnomAD, the gene has a low tolerance for loss-of-function variants (pLI = 1, o/e = 0.07). We functionally evaluated the PPP5C variant in C. elegans by knocking the variant into the orthologous gene, pph-5, at the corresponding residue, Ala48Thr. We employed assays in three different biological processes where pph-5 was known to function through opposing the activity of genes, mec-15 and sep-1. We demonstrated that, in contrast to control animals, the pph-5 Ala48Thr variant suppresses the neurite growth phenotype and the GABA signaling defects of mec-15 mutants, and the embryonic lethality of sep-1 mutants. The Ala48Thr variant did not display dominance and behaved similarly to the reference pph-5 null, indicating that the variant is likely a strong hypomorph or complete loss-of-function. We conclude that pph-5 Ala48Thr is damaging in C. elegans. By extension in the proband, PPP5C p.Ala47Thr is likely damaging, the de novo dominant presentation is consistent with haplo-insufficiency, and the PPP5C variant is likely responsible for one or more of the proband's phenotypes.

A dominant negative variant of RAB5B disrupts maturation of surfactant protein B and surfactant protein C.

Pathogenic variants in surfactant proteins SP-B and SP-C cause surfactant deficiency and interstitial lung disease. Surfactant proteins are synthesized as precursors (proSP-B, proSP-C), trafficked, and processed via a vesicular-regulated secretion pathway; however, control of vesicular trafficking events is not fully understood. Through the Undiagnosed Diseases Network, we evaluated a child with interstitial lung disease suggestive of surfactant deficiency. Variants in known surfactant dysfunction disorder genes were not found in trio exome sequencing. Instead, a de novo heterozygous variant in RAB5B was identified in the Ras/Rab GTPases family nucleotide binding domain, p.Asp136His. Functional studies were performed in Caenorhabditis elegans by knocking the proband variant into the conserved position (Asp135) of the ortholog, rab-5 Genetic analysis demonstrated that rab-5[Asp135His] is damaging, producing a strong dominant negative gene product. rab-5[Asp135His] heterozygotes were also defective in endocytosis and early endosome (EE) fusion. Immunostaining studies of the proband's lung biopsy revealed that RAB5B and EE marker EEA1 were significantly reduced in alveolar type II cells and that mature SP-B and SP-C were significantly reduced, while proSP-B and proSP-C were normal. Furthermore, staining normal lung showed colocalization of RAB5B and EEA1 with proSP-B and proSP-C. These findings indicate that dominant negative-acting RAB5B Asp136His and EE dysfunction cause a defect in processing/trafficking to produce mature SP-B and SP-C, resulting in interstitial lung disease, and that RAB5B and EEs normally function in the surfactant secretion pathway. Together, the data suggest a noncanonical function for RAB5B and identify RAB5B p.Asp136His as a genetic mechanism for a surfactant dysfunction disorder.

Sex- and mutation-specific p53 gain-of-function activity in gliomagenesis.

In cancer, missense mutations in the DNA-binding domain of TP53 are common. They abrogate canonical p53 activity and frequently confer gain-of-oncogenic function (GOF) through localization of transcriptionally active mutant p53 to non-canonical genes. We found that several recurring p53 mutations exhibit a sex difference in frequency in patients with glioblastoma (GBM). In vitro and in vivo analysis of three mutations, p53R172H, p53Y202C, and p53Y217C revealed unique interactions between cellular sex and p53 GOF mutations that determined each mutation's ability to transform male versus female primary mouse astrocytes. These phenotypic differences were correlated with sex- and p53 mutation- specific patterns of genomic localization to the transcriptional start sites of upregulated genes belonging to core cancer pathways. The promoter regions of these genes exhibited a sex difference in enrichment for different transcription factor DNA-binding motifs. Together, our data establish a novel mechanism for sex specific mutant p53 GOF activity in GBM with implications for all cancer.

Effect of Whole-Genome Sequencing on the Clinical Management of Acutely Ill Infants With Suspected Genetic Disease: A Randomized Clinical Trial.

Importance: Whole-genome sequencing (WGS) shows promise as a first-line genetic test for acutely ill infants, but widespread adoption and implementation requires evidence of an effect on clinical management. Objective: To determine the effect of WGS on clinical management in a racially and ethnically diverse and geographically distributed population of acutely ill infants in the US. Design, Setting, and Participants: This randomized, time-delayed clinical trial enrolled participants from September 11, 2017, to April 30, 2019, with an observation period extending to July 2, 2019. The study was conducted at 5 US academic medical centers and affiliated children's hospitals. Participants included infants aged between 0 and 120 days who were admitted to an intensive care unit with a suspected genetic disease. Data were analyzed from January 14 to August 20, 2020. Interventions: Patients were randomized to receive clinical WGS results 15 days (early) or 60 days (delayed) after enrollment, with the observation period extending to 90 days. Usual care was continued throughout the study. Main Outcomes and Measures: The main outcome was the difference in the proportion of infants in the early and delayed groups who received a change of management (COM) 60 days after enrollment. Additional outcome measures included WGS diagnostic efficacy, within-group COM at 90 days, length of hospital stay, and mortality. Results: A total of 354 infants were randomized to the early (n = 176) or delayed (n = 178) arms. The mean participant age was 15 days (IQR, 7-32 days); 201 participants (56.8%) were boys; 19 (5.4%) were Asian; 47 (13.3%) were Black; 250 (70.6%) were White; and 38 (10.7%) were of other race. At 60 days, twice as many infants in the early group vs the delayed group received a COM (34 of 161 [21.1%; 95% CI, 15.1%-28.2%] vs 17 of 165 [10.3%; 95% CI, 6.1%-16.0%]; P = .009; odds ratio, 2.3; 95% CI, 1.22-4.32) and a molecular diagnosis (55 of 176 [31.0%; 95% CI, 24.5%-38.7%] vs 27 of 178 [15.0%; 95% CI, 10.2%-21.3%]; P < .001). At 90 days, the delayed group showed a doubling of COM (to 45 of 161 [28.0%; 95% CI, 21.2%-35.6%]) and diagnostic efficacy (to 56 of 178 [31.0%; 95% CI, 24.7%-38.8%]). The most frequent COMs across the observation window were subspecialty referrals (39 of 354; 11%), surgery or other invasive procedures (17 of 354; 4%), condition-specific medications (9 of 354; 2%), or other supportive alterations in medication (12 of 354; 3%). No differences in length of stay or survival were observed. Conclusions and Relevance: In this randomized clinical trial, for acutely ill infants in an intensive care unit, introduction of WGS was associated with a significant increase in focused clinical management compared with usual care. Access to first-line WGS may reduce health care disparities by enabling diagnostic equity. These data support WGS adoption and implementation in this population. Trail Registration: ClinicalTrials.gov Identifier: NCT03290469.

Functional analysis of a de novo variant in the neurodevelopment and generalized epilepsy disease gene NBEA.

Neurobeachin (NBEA) was initially identified as a candidate gene for autism. Recently, variants in NBEA have been associated with neurodevelopmental delay and childhood epilepsy. Here, we report on a novel NBEA missense variant (c.5899G > A, p.Gly1967Arg) in the Domain of Unknown Function 1088 (DUF1088) identified in a child enrolled in the Undiagnosed Diseases Network (UDN), who presented with neurodevelopmental delay and seizures. Modeling of this variant in the Caenorhabditis elegans NBEA ortholog, sel-2, indicated that the variant was damaging to in vivo function as evidenced by altered cell fate determination and trafficking of potassium channels in neurons. The variant effect was indistinguishable from that of the reference null mutation suggesting that the variant is a strong hypomorph or a complete loss-of-function. Our experimental data provide strong support for the molecular diagnosis and pathogenicity of the NBEA p.Gly1967Arg variant and the importance of the DUF1088 for NBEA function.

Altered neuronal physiology, development, and function associated with a common chromosome 15 duplication involving CHRNA7.

BACKGROUND: Copy number variants (CNVs) linked to genes involved in nervous system development or function are often associated with neuropsychiatric disease. While CNVs involving deletions generally cause severe and highly penetrant patient phenotypes, CNVs leading to duplications tend instead to exhibit widely variable and less penetrant phenotypic expressivity among affected individuals. CNVs located on chromosome 15q13.3 affecting the alpha-7 nicotinic acetylcholine receptor subunit (CHRNA7) gene contribute to multiple neuropsychiatric disorders with highly variable penetrance. However, the basis of such differential penetrance remains uncharacterized. Here, we generated induced pluripotent stem cell (iPSC) models from first-degree relatives with a 15q13.3 duplication and analyzed their cellular phenotypes to uncover a basis for the dissimilar phenotypic expressivity. RESULTS: The first-degree relatives studied included a boy with autism and emotional dysregulation (the affected proband-AP) and his clinically unaffected mother (UM), with comparison to unrelated control models lacking this duplication. Potential contributors to neuropsychiatric impairment were modeled in iPSC-derived cortical excitatory and inhibitory neurons. The AP-derived model uniquely exhibited disruptions of cellular physiology and neurodevelopment not observed in either the UM or unrelated controls. These included enhanced neural progenitor proliferation but impaired neuronal differentiation, maturation, and migration, and increased endoplasmic reticulum (ER) stress. Both the neuronal migration deficit and elevated ER stress could be selectively rescued by different pharmacologic agents. Neuronal gene expression was also dysregulated in the AP, including reduced expression of genes related to behavior, psychological disorders, neuritogenesis, neuronal migration, and Wnt, axonal guidance, and GABA receptor signaling. The UM model instead exhibited upregulated expression of genes in many of these same pathways, suggesting that molecular compensation could have contributed to the lack of neurodevelopmental phenotypes in this model. However, both AP- and UM-derived neurons exhibited shared alterations of neuronal function, including increased action potential firing and elevated cholinergic activity, consistent with increased homomeric CHRNA7 channel activity. CONCLUSIONS: These data define both diagnosis-associated cellular phenotypes and shared functional anomalies related to CHRNA7 duplication that may contribute to variable phenotypic penetrance in individuals with 15q13.3 duplication. The capacity for pharmacological agents to rescue some neurodevelopmental anomalies associated with diagnosis suggests avenues for intervention for carriers of this duplication and other CNVs that cause related disorders.

Model organisms contribute to diagnosis and discovery in the undiagnosed diseases network: current state and a future vision.

Decreased sequencing costs have led to an explosion of genetic and genomic data. These data have revealed thousands of candidate human disease variants. Establishing which variants cause phenotypes and diseases, however, has remained challenging. Significant progress has been made, including advances by the National Institutes of Health (NIH)-funded Undiagnosed Diseases Network (UDN). However, 6000-13,000 additional disease genes remain to be identified. The continued discovery of rare diseases and their genetic underpinnings provides benefits to affected patients, of whom there are more than 400 million worldwide, and also advances understanding the mechanisms of more common diseases. Platforms employing model organisms enable discovery of novel gene-disease relationships, help establish variant pathogenicity, and often lead to the exploration of underlying mechanisms of pathophysiology that suggest new therapies. The Model Organism Screening Center (MOSC) of the UDN is a unique resource dedicated to utilizing informatics and functional studies in model organisms, including worm (Caenorhabditis elegans), fly (Drosophila melanogaster), and zebrafish (Danio rerio), to aid in diagnosis. The MOSC has directly contributed to the diagnosis of challenging cases, including multiple patients with complex, multi-organ phenotypes. In addition, the MOSC provides a framework for how basic scientists and clinicians can collaborate to drive diagnoses. Customized experimental plans take into account patient presentations, specific genes and variant(s), and appropriateness of each model organism for analysis. The MOSC also generates bioinformatic and experimental tools and reagents for the wider scientific community. Two elements of the MOSC that have been instrumental in its success are (1) multidisciplinary teams with expertise in variant bioinformatics and in human and model organism genetics, and (2) mechanisms for ongoing communication with clinical teams. Here we provide a position statement regarding the central role of model organisms for continued discovery of disease genes, and we advocate for the continuation and expansion of MOSC-type research entities as a Model Organisms Network (MON) to be funded through grant applications submitted to the NIH, family groups focused on specific rare diseases, other philanthropic organizations, industry partnerships, and other sources of support.

Biallelic ASCC1 variants including a novel intronic variant result in expanded phenotypic spectrum of spinal muscular atrophy with congenital bone fractures 2 (SMABF2).

Spinal muscular atrophy with congenital bone fractures 2 (SMABF2), a type of arthrogryposis multiplex congenita (AMC), is characterized by congenital joint contractures, prenatal fractures of long bones, and respiratory distress and results from biallelic variants in ASCC1. Here, we describe an infant with severe, diffuse hypotonia, congenital contractures, and pulmonary hypoplasia characteristic of SMABF2, with the unique features of cleft palate, small spleen, transverse liver, and pulmonary thromboemboli with chondroid appearance. This infant also had impaired coagulation with diffuse petechiae and ecchymoses which has only been reported in one other infant with AMC. Using trio whole genome sequencing, our proband was identified to have biallelic variants in ASCC1. Using deep next generation sequencing of parental cDNA, we characterized alteration of splicing encoded by the novel, maternally inherited ASCC1 variant (c.297-8 T > G) which provides a mechanism for functional pathogenicity. The paternally inherited ASCC1 variant is a rare nonsense variant (c.466C > T; p.Arg156*) that has been previously identified in one other infant with AMC. This report extends the phenotypic characteristics of ASCC1-associated AMC (SMABF2) and describes a novel intronic variant that partially disrupts RNA splicing.

Commonalities across computational workflows for uncovering explanatory variants in undiagnosed cases.

PURPOSE: Genomic sequencing has become an increasingly powerful and relevant tool to be leveraged for the discovery of genetic aberrations underlying rare, Mendelian conditions. Although the computational tools incorporated into diagnostic workflows for this task are continually evolving and improving, we nevertheless sought to investigate commonalities across sequencing processing workflows to reveal consensus and standard practice tools and highlight exploratory analyses where technical and theoretical method improvements would be most impactful. METHODS: We collected details regarding the computational approaches used by a genetic testing laboratory and 11 clinical research sites in the United States participating in the Undiagnosed Diseases Network via meetings with bioinformaticians, online survey forms, and analyses of internal protocols. RESULTS: We found that tools for processing genomic sequencing data can be grouped into four distinct categories. Whereas well-established practices exist for initial variant calling and quality control steps, there is substantial divergence across sites in later stages for variant prioritization and multimodal data integration, demonstrating a diversity of approaches for solving the most mysterious undiagnosed cases. CONCLUSION: The largest differences across diagnostic workflows suggest that advances in structural variant detection, noncoding variant interpretation, and integration of additional biomedical data may be especially promising for solving chronically undiagnosed cases.

Discovery of a novel CHD7 CHARGE syndrome variant by integrated omics analyses.

Chromodomain helicase DNA-binding protein 7 (CHD7) pathogenic variants are identified in more than 90% of infants and children with CHARGE (Coloboma of the iris, retina, and/or optic disk; congenital Heart defects, choanal Atresia, Retardation of growth and development, Genital hypoplasia, and characteristic outer and inner Ear anomalies and deafness) syndrome. Approximately, 10% of cases have no known genetic cause identified. We report a male child with clinical features of CHARGE syndrome and nondiagnostic genetic testing that included chromosomal microarray, CHD7 sequencing and deletion/duplication analysis, SEMA3E sequencing, and trio exome and whole-genome sequencing (WGS). We used a comprehensive clinical assessment, genome-wide methylation analysis (GMA), reanalysis of WGS data, and CHD7 RNA studies to discover a novel variant that causes CHD7 haploinsufficiency. The 7-year-old Hispanic male proband has typical phenotypic features of CHARGE syndrome. GMA revealed a CHD7-associated epigenetic signature. Reanalysis of the WGS data with focused bioinformatic analysis of CHD7 detected a novel, de novo 15 base pair deletion in Intron 4 of CHD7 (c.2239-20_2239-6delGTCTTGGGTTTTTGT [NM_017780.3]). Using proband RNA, we confirmed that this novel deletion causes CHD7 haploinsufficiency by disrupting the canonical 3' splice site and introducing a premature stop codon. Integrated genomic, epigenomic, and transcriptome analyses discovered a novel CHD7 variant that causes CHARGE syndrome.

Precise breakpoint detection in a patient with 9p- syndrome.

We present a case of 9p- syndrome with a complex chromosomal event originally characterized by the classical karyotype approach as 46,XX,der(9)t(9;13)(p23;q13). We used advanced technologies (Bionano Genomics genome imaging and 10× Genomics sequencing) to characterize the location of the translocation and accompanying deletion on Chromosome 9 and duplication on Chromosome 13 with single-nucleotide breakpoint resolution. The translocation breakpoint was at Chr 9:190938 and Chr 13:50850492, the deletion at Chr 9:1-190938, and the duplication at Chr 13:50850492-114364328. We identified genes in the deletion and duplication regions that are known to be associated with this patient's phenotype (e.g., ZIC2 in dysmorphic facial features, FOXD4 in developmental delay, RNASEH2B in developmental delay, and PCDH9 in autism). Our results indicate that clinical genomic assessment of individuals with complex karyotypes can be refined to a single-base-pair resolution when utilizing Bionano and 10× Genomics sequencing. With the 10× Genomics data, we were also able to characterize other variation (e.g., loss of function) throughout the remainder of the patient's genome. Overall, the Bionano and 10× technologies complemented each other and provided important insight into our patient with 9p- syndrome. Altogether, these results indicate that newer technologies can identify precise genomic variants associated with unique patient phenotypes that permit discovery of novel genotype-phenotype correlations and therapeutic strategies.

MYT1L: A systematic review of genetic variation encompassing schizophrenia and autism.

Variations in MYT1L, a gene encoding a transcription factor expressed in the brain, have been associated with autism, intellectual disability, and schizophrenia. Here we provide an updated review of published reports of neuropsychiatric correlates of loss of function and duplication of MYT1L. Of 27 duplications all were partial; 33% were associated exclusively with schizophrenia, and the chromosomal locations of schizophrenia-associated duplications exhibited a distinct difference in pattern-of-location from those associated with autism and/or intellectual disability. Of 51 published heterozygous loss of function variants, all but one were associated with intellectual disability, autism, or both, and one resulted in no neuropsychiatric diagnosis. There were no reports of schizophrenia associated with loss of function variants of MYT1L (Fisher's exact p < .00001, for contrast with all reported duplications). Although the precise function of the various mutations remains unspecified, these data collectively establish the candidacy of MYT1L as a reciprocal mutation, in which schizophrenia may be engendered by partial duplications, typically involving the 3' end of the gene, while developmental disability-notably autism-is associated with both loss of function and partial duplication. Future research on the specific effects of contrasting mutations in MYT1L may provide insight into the causal origins of autism and schizophrenia.

Phenotypic expansion of KMT2D-related disorder: Beyond Kabuki syndrome.

Pathogenic variants in KMT2D, which encodes lysine specific methyltransferase 2D, cause autosomal dominant Kabuki syndrome, associated with distinctive dysmorphic features including arched eyebrows, long palpebral fissures with eversion of the lower lid, large protuberant ears, and fetal finger pads. Most disease-causing variants identified to date are putative loss-of-function alleles, although 15-20% of cases are attributed to missense variants. We describe here four patients (including one previously published patient) with de novo KMT2D missense variants and with shared but unusual clinical findings not typically seen in Kabuki syndrome, including athelia (absent nipples), choanal atresia, hypoparathyroidism, delayed or absent pubertal development, and extreme short stature. These individuals also lack the typical dysmorphic facial features found in Kabuki syndrome. Two of the four patients had severe interstitial lung disease. All of these variants cluster within a 40-amino-acid region of the protein that is located just N-terminal of an annotated coiled coil domain. These findings significantly expand the phenotypic spectrum of features associated with variants in KMT2D beyond those seen in Kabuki syndrome and suggest a possible new underlying disease mechanism for these patients.